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1 Introduction 



The collective transverse in-plane flow has been used extensively over the past three 
decades to study the properties of hot and dense nuclear matter., i.e., the nuclear matter 
equation of state (EOS) and in- medium nucleon-nucleon cross section [1]. It has been 
reported to be highly sensitive to the above mentioned properties and also to the entrance 
channel parameters like incident energy, colliding geometry and system size [2-6]. The 
energy dependence of flow led to its disappearance at a particular incident energy called 
energy of vanishing flow (EVP) or balance energy {Ef,ai) [7]. A large number of theoretical 
studies have been carried out in the past studying the sensitivity of E^^^ to the system 
size and coUiding geometry [8-10]. 

Role of isospin degree of freedom in collective transverse in-plane flow and its dis- 
appearance has also been a matter of great interest for the past decade [11, 12]. Isospin 
degree of freedom plays its role in determining the nuclear equation of state of asymmetric 
nuclear matter. The availability of radioactive ion beams (RIBs) [13, 14] around the world 
helps in carrying out the studies on the matter lying far away from the stability line. A 
number of studies have been carried out in the recent past to see the role of isospin degree 
of freedom in collective flow and its disappearance [11, 12, 15]. In Ref. [16] author and 
others studied the isospin effects in Ei,ai at all the coUiding geometries. A very few studies 
have been carried out to study other related phenomena at E^ai of the neutron-rich sys- 
tems. In Ref. [17] Sood and Puri have presented a complete study of the nuclear dynamics 
at Ehai for stable systems. The study includes participant-spectator matter, density and 
temperature reached in a heavy-ion reaction at E^ai- Motivated by this, author and others 
presented a study of participant-spectator matter of neutron-rich systems at E^a/ in Ref. 
[18]. The study revealed a similar behaviour of participant-spectator matter for neutron- 
rich systems as for the stable systems. Another important quantity which reflects the 
dynamics in a heavy-ion collision is the density and temperature reached in a reaction. In 
the present paper, we study the density and temperature reached in heavy-ion reactions 
of neutron- rich matter at Efjai- We also aim to see the role of isospin degree of freedom 
in the density and temperature reached in the reactions of neutron-rich systems and to 
see if the behaviour of these properties at balance energy differs from that for systems lie 
close to the stability line. 

The present study is carried out within the framework of isospin- dependent quantum 
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molecular dynamics (IQMD) model [19]. Section 2 describes the model in brief. Section 
3 explains the results and gives our discussion, and Sec. 4 summarizes the results. 

2 The model 

The IQMD model [19] which is the extension of quantum molecular dynamics (QMD) [20] 
model treats different charge states of nucleons, deltas, and pions explicitly, as inherited 
from the Vlasov-Uehling-Uhlenbeck (VUU) model. The IQMD model has been used suc- 
cessfully for the analysis of a large number of observables from low to relativistic energies. 
Puri and coworkers have demonstrated that QMD, IQMD carries essential physics needed 
to demonstrate the various phenomena such as collective flow, multifragmentation and 
particle production [21, 22]. The isospin degree of freedom enters into the calculations via 
symmetry potential, cross sections, and Coulomb interaction. 

In this model, baryons are represented by Gaussian-shaped density distributions 

1 1 2L 

= ;;^^''P^~[''~'^^^)]'^) ^ exp(-[p-pi(t)]2— ) (1) 

Nucleons are initialized in a sphere with radius R = 1.12 A^/^ fm, in accordance with 
liquid-drop model. Each nucleon occupies a volume of /i^, so that phase space is uniformly 
filled. The initial momenta are randomly chosen between and Fermi momentum {pf)- 
The nucleons of the target and projectile interact by two- and three-body Skyrme forces, 
Yukawa potential, and Coulomb interactions. In addition to the use of explicit charge 
states of all baryons and mesons, a symmetry potential between protons and neutrons 
corresponding to the Bethe-Weizsacker mass formula has been included. The hadrons 
propagate using Hamilton equations of motion: 

dfj _ d{H) _ dpi _ d{H) 
dt dpi ' dt dfi 

with 



(H) = {T) + {V) 



i j>i 



Mf,p,t)V'^if',r) 



xfj{r',p',t)dfdf' 



(3) 
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The baryon potential V^j in the above relation, reads as 



V^'j(f'-f) 



Skyrmc 



Yukawa 




sym 



[M(r '-f)+M(r '-r)p^-^( 



)] 



2 



exp(|(f'-r)|///) , Z,Z,-e2 



(|(f'-r)|//.) |(f'-r)| 



+i4-T3iT3j5(rl'-r;). 



(4) 



Here and Zj denote the charges of ith and jth baryon, and T^i and are their 
respective Ts components (i.e., 1/2 for protons and —1/2 for neutrons). The parameters ji 
and ti,...,ti are adjusted to the real part of the nucleonic optical potential. For the density 
dependence of the nucleon optical potential, standard Skyrme type parametrization is 
employed. Wc use a soft equation of state along with the standard isospin- and energy- 
dependent cross section reduced by 20%, i.e. a = 0.8 crl^^^"^. In a recent study, Gautam et 
al. [15] has confronted the theoretical calculations of IQMD with the data of ^^Ni Ni 
and ^^Fe+^^Fe [11]. The results with the soft EOS (along with the momentum- dependent 
interactions) and above choice of cross section are in good agreement with the data at 
all colliding geometries. The details about the elastic and inelastic cross sections for 
proton-proton and proton-neutron collisions can be found in [19,23]. The cross sections 
for neutron-neutron collisions are assumed to be equal to the proton-proton cross sections. 
Two particles collide if their minimum distance d fulfills 



where 'type' denotes the ingoing coUision partners (N-N....). Exphcit Pauh blocking is 
also included; i.e. Pauh blocking of the neutrons and protons is treated separately. We 
assume that each nucleon occupies a sphere in coordinate and momentum space. This 
trick yields the same Pauh blocking ratio as an exact calculation of the overlap of the 
Gaussians will yield. We calculate the fractions Pi and P2 of final phase space for each 
of the two scattering partners that arc already occupied by other nucleons with the same 
isospin as that of scattered ones. The collision is blocked with the probability 




(5) 



Puock = 1 - [1 - min{Pi, 1)][1 - min{P2, 1)], 



(6) 
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and, correspondingly is allowed with the probability 1 - Fuock- For a nucleus in its 
ground state, we obtain an averaged blocking probability (Puock) = 0.96. Whenever an 
attempted collision is blocked, the scattering partners maintain the original momenta 
prior to scattering. 

3 Results and discussion 

We simulate the reactions of Ca+Ca, Ni+Ni, Zr+Zr, Sn+Sn, and Xe+Xe series hav- 
ing N/Z = 1.0, 1.6 and 2.0. In particular, we simulate the reactions of ^'^Ca+^'^Ca 
(105), 52Ca+52c^ (g5)^ 60ca+60ca (73); ^^Ni+ssNi (98), ^^Ni+^^^j (32), 84^1+84^1 (72); 
s^Zr+s^Zr (86), ^o^Zr+^o^Zr (74), ^^ozr+i^ozr (67); i°0Sn+i°0Sn (82), '^^Sn+'^^Sn (72), 
i50gn^i50gn (64) and iioXe+iioXe (76), ^^oxe+i^oxe (68) and ^^^Xe+^^^Xe (61) at an im- 
pact parameter of b/bmax = 0.2-0.4 and at the incident energies equal to balance energy. 
The values in the brackets represent the balance energies for the systems. The reactions 
are followed till the transverse in-plane flow saturates. It is worth mentioning here that 
the saturation time varies with the mass of the system. It has been shown in Ref. [24] 
that the transverse in-plane flow in lighter colliding nuclei saturates earlier compared to 
heavy coUiding nuclei. Saturation time is about 100 (150 fm/c) in hghter (heavy) colliding 
nuclei in the present energy domain. We use the quantity " directed transverse momentum 
{pf^)" to define the nuclear transverse in-plane flow, which is deflned as [19, 20, 25] 

1 

^T) = JT.''9n{ymPo.{^), (7) 

i=l 

where y{i) andpa;(i) are, respectively, the rapidity (calculated in the center of mass system) 
and the momentum of the i*^ particle. The rapidity is deflned as 

y(,) = i,n|W±M), (8) 
2 E(i)-pji)' 

where E{i) and Pzii) are, respectively, the energy and longitudinal momentum of the 
i*^ particle. In this definition, all the rapidity bins arc taken into account. 

In fig. 1(a), wc display the time evolution of average density (p"'"^/ po) whereas fig. 1(b) 
displays the time evolution of maximum density {p"^""^ / po) for the systems having N/Z 
= 1.0, i.e, we display the reactions of ^oCa+^oCa, ^^Ni+^^Ni, ^^Zr+^^Zr, looSn+^ooSn, 
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Time (fm/c) 

Figure 1: (Color online) The time evolution of spectator matter (left panels) and partici- 
pant matter (right panels) for systems having N/Z = 1.0, 1.6 and 2.0. Lines are explained 
in the text. 

and ^^^Xe+^^^Xe at energy equal to balance energy. Lines represent the different sys- 
tems. Solid, dashed, dotted, dash-dotted, and short dotted lines represent the reactions 
of Ca+Ca, Ni+Ni, Zr+Zr, Sn+Sn, and Xe-|-Xe, respectively. From figure, we find that 
maximal value of p""^^ / po is higher for lighter systems as compared to the heavier ones. 
Moreover, the density profile is more extended in heavier systems indicating that the 
reaction finishes later in heavier systems. This is because of the fact that the heavier 
reaction occurs at low incident energy. Also the p""^^ / Po and p^""^ / po are nearly same 
for heavier systems but differ for lighter systems as in Ref. Further, the maximum and 
average densities are comparable for medium and heavy mass systems indicating that the 
dense matter is formed widely and uniformly in the central zone of the reaction. On the 
other hand, the substantial difference in two densities for the lighter colliding nuclei has 
been explained in Ref. [17] and indicates the non-homogeneous nature of dense matter. It 
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Figure 2: (Color online) The N/Z dependence of participant and spectator matter. Sym- 
bols are explained in the text. 



is worth mentioning that collective flow saturates at higher densities whereas multifrag- 
mentation occurs at sub-density zone. Other phenomena such as fusion-fission are very 
low density phenomena [26]. 

The quantity which reflects the density achieved in a reaction is the collision rate. In 
flg. 2, we display the time evolution of the collision rate for various systems having N/Z = 
1.0, 1.6 and 2.0. Solid, dashed and dotted lines corresponds to systems having N/Z = 1.0, 
1.6 and 2.0, respectively. From flgure, we see that collision rate flrst increases with time, 
reaches maximum at around 20-40 fm/c (which is the high dense phase of the reaction) 
and then flnally decreases and becomes constant at around 80 fm/c. We also flnd that 
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Figure 3: The system size dependence of participant and spectator matter for different 
N/Z ratios. Various symbols are explained in the text. 

the maximum value of the collision rate also increases with the system mass. Moreover, 
the effect of N/Z ratio on the coUision rate is very less. 

In fig. 3 we display the system size dependence of maximal value of the maximum 
^pinax^ and average density {p"-^^) for the systems having N/Z = 1.0, 1.6 and 2.0. We see 
that the maximal value of p^""^ and p"^^ follows a power law behaviour proportional to 



8 





2.0 




1.8 




1.6 




1.4 




1.2 




1.0 




2.2 




2.0 






Q. 


1.8 


A 


1.6 




1.4 


V 


1.2 




1.0 




2.2 




2.0 




1.8 




1.6 




1.4 




1.2 




1.0 



'Ml 

- O 


o 


I'M 
O - 


'_ • 


• 


• _ 


M , 


1 , 1 


Ca+Ca! 

Ml, 



MM 


1 ' M M M 


- O 




~ • 






• - 




Ni+Ni: 




'Ml 


1 M ' M M 


- O 


O 1 


r • 


• o - 






"(c) 


Zr+Zr: 




MM 


1 ' M M M 


-_ O 


O 

• o : 
• - 


- • 




Sn+Sn: 




2.2 
2.0 
1.8 
1.6 
1.4 
1.2 
1.0 

2.2 
2.0 
1.8 
1.6 
1.4 
1.2 
1.0 



'Ml 


1 M ' M M 




'_ O 






-_ • 








O <P>"" 


:(e) 

.1,1 


Xe+Xe" 
1,1,1,1, 





0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 

N/Z 



Figure 4: The N/Z dependence of participant and spectator matter. Symbols are ex- 
plained in the text. 

A^. The power law factor is O.Oli 0.04 (-0.03± 0.04), -0.04± 0.01 (-0.08± 0.02), and 
-0.07± 0.008 (-0.10± 0.02) for p""'f ( p™"^) for systems with N/Z = 1.0, 1.6 and 2.0, 
respectively. It shows that the dependence of maximal value of p""""^ and p™-"^ is very 
weak on the system size for all the N/Z ratios. This was also predicted in Ref.citesood2 
for stable systems. 

In fig. 4 we display the N/Z dependence of maximal value of p"''"^ and p^"'^. Sohd 
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Figure 5: The system size dependence of anisotropy ratio for various N/Z ratios. 

(open) symbols display the results for p""^ ( p"^"^). Prom figure we see that both p""^ 
and p"''"^ decreases shghtly with N/Z of the system for all the system masses. A slight 
exception to this is there for the lighter mass of Ca+Ca. 

The another associated quantity linked with the dense matter is the temperature. 
In principle, a true temperature can be defined only for a thermalized and equilibrated 
matter. Since in heavy-ion collisions the matter is non-equilibrated, one can not talk of 
"temperature". One can, however, look in terms of the local environment only. In our 
present case, we follow the description of the temperature given in Refs. [27, 28]. In the 
present case, extraction of the temperature Tis based on the local density approximation. 
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Figure 6: The system size dependence of anisotropy ratio for various N/Z ratios. 

i.e., one deduces the temperature in a volume element surrounding the position of each 
particle at a given time step [27, 28] . Here, we postulate that each local volume element 
of nuclear matter in coordinate space and time has some "temperature" defined by the 
diffused edge of the deformed Fermi distribution consisting of two coUiding Fermi spheres, 
which is typical for a nonequilibrium momentum distribution in heavy-ion collisions. 

In this formahsm (dubbed the hot Thomas- Fermi approach [27]), one determines ex- 
tensive quantities like the density and kinetic energy as well as entropy with the help of 
momentum distributions at a given temperature. Using this formalism, we also extracted 
the average and maximum temperature within a central sphere of 2 fm radius as described 
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in the case of density. 

In fig. 5 we display the maximal value of < T""^ > and < T™""^ > as a function of the 
composite mass of the system. From figure, we see that < T"^^ > and < T"*"^ > follows 
a power law behaviour with system mass for all the N/Z ratios. The power law factor 
is -0.16± 0.06 (-0.08± 0.05), -0.15± 0.07 (-0.16± 0.06), and -0.19± 0.09 (-0.15± 0.07) 
for for < > (< T"^"^ >) for systems having N/Z = 1.0, 1.6 and 2.0, respectively. 
Similar power law behaviour was also predicted in Ref. [17] for stable systems (N/Z ~ 1). 
This system size dependence of temperature is sharp in contrast with the density. This 
is because of the fact that the temperature depends on the kinetic energy of the system. 

In fig. 6 we display the N/Z dependence of maximal value of < T"-"^ > and < T"^"-^ > 
for various system masses. Solid (open) symbols represent the maximal value of < T"-'"^ > 
(< T*""^ >). From figure, we see that for all the system masses, < T"''"^ > and < T™"^ > 
decreases with N/Z of the system. 

4 Summary 

We studied the collision rate, density and temperature reached in reactions of neutron- 
rich systems at energy of vanishing flow. Our results pointed the similar behvaiour for 
neutron-rich systems as for the stable systems. We also investigated the mass dependence 
of these quantities. We found a very weak mass dependence of the density although the 
temperature follows a significant mass dependence. 

This work has been supported by a grant from Centre of Scientific and Industrial 
Research (CSIR), Govt, of India. Author is thankful to Profs. J. Aichelin and R. K. Puri 
for enlightening discussions on the present work. 
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